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ABSTRACT 

In  this  paper  we  analyze  the  behavior  of  the  so-called  p-version  of  the 
finite  element  method  when  applied  to  the  equations  of  plane  strain  linear 
elasticity.  We  establish  optimal  rate  error  estimates  that  are  uniformly 
valid,  independent  of  the  value  of  the  Poisson  ratio,  v,  in  the  interval 
]0,1/2[.  This  shows  that  the  p-version  does  not  exhibit  the  degeneracy 
phenomenon  which  has  led  to  the  use  of  various,  only  partially  justified 
techniques  of  reduced  integration  or  mixed  formulations  for  more  standard 
finite  element  schemes  and  the  case  of  a  nearly  incompressible  material# 
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1 .  Introduction 

It  is  well  known  that  standard  finite  element  schemes 
based  on  elements  and  direct  minimization  of  the  displacement 
energy  have  very  poor  convergence  properties  for  nearly 
incompressible  materials.  The  same  degeneracy  occurs  in 
connection  with  plate-bending  models  that  account  for  shear, 
such  as  the  Timoshenko  plate  (or  beam)  model,  when  the 
thickness  of  the  structure  is  small  (cf.  [17]).  In  order 
to  get  reliable  results  one  would  have  to  take  a  mesh  size  h 
which  in  the  first  case  is  much  less  than  1-2V  (v  denotes  the 
Poisson  ratio)  and  in  the  second  case  much  less  than  the  thick¬ 
ness  of  the  structure.  This  is  in  most  cases  not  practically 
feasible. 

One  way  to  try  to  avoid  the  latter  problem  is  to  include 
more  regular  elements,  e.g.  (cf .  [11] ) ,  but  in  more  than  one 

dimension  the  extra  continuity  conditions  are  very  difficult  to 
handle  in  practice. 

By  far  the  most  simple  and  successful  remedy  in  both  cases 
has  been  the  technique  of  "reduced  integration"  (cf.  [18]).  The 
idea  is  to  replace  some  of  the  integrals  in  the  energy  by  the 
corresponding  expressions  arising  from  the  application  of  some 
quadrature  rule.  The  finite  element  spaces  remain  simple  C^,  for 
example,  piecewise  linear.  This  approach  of  course  is  only  a  partial 
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remedy  since  it  cannot  provide  better  results  measured  in  terms 
of  energy,  but  it  makes  it  possible  to  obtain  good  approximations 
to  the  individual  components  of  the  displacement.  In  many 
instances  there  is  a  close  connection  between  this  approach 
of  reduced  integration  and  a  so-called  mixed  formulation  of 
the  variational  problem  (cf.  [2,14]).  There  are  numerous 
papers  in  the  engineering  literature  devoted  to  the  questions 
of  what  terms  in  the  energy  to  evaluate  only  approximately 
and  what  quadrature  rule  to  use  for  a  particular  element. 

It  is  quite  clear  however  that  no  comprehensive  mathematical 
analysis  exists  at  this  time  to  justify  the  benefits  of 
reduced  integration.  Some  partial  justification  may  be  found 
e.g.  in  [2] ,  for  the  case  of  the  Timoshenko  beam  model  and  small 
thickness,  and  in  [12]  for  the  case  of  one  particular  element  and 
either  an  almost  incompressible  material  or  the  Timoshenko 
plate  model. 

Finite  element  methods  in  general  are  based  on  some 
partition  of  the  domain  (typically  a  triangulation)  and  function 
spaces  whose  elements  are  piecewise  polynomials  on  each  triangle. 
The  h-version  of  the  finite  element  method  (the  one  discussed 
so  far)  obtains  convergent  solutions  by  letting  the  mesh  size  h 
(i.e.,  the  diameter  of  the  largest  triangle)  tend  to  0  using 
only  piecewise  polynomials  of  some  fixed  degree.  In  contrast 
the  so-called  p-version  of  the  finite  element  method  is  based  on 
a  fixed  triangulation  and  the  degree  of  the  piecewise  polynomials 
is  increased  in  order  to  give  convergence.  For  more  details  about 
the  convergence  properties  of  the  p-version  and  combinations  of 
the  two  versions  see  [5,6,7]. 


In  this  paper  we  study  the  behavior  of  the  p-version  for 
an  almost  incompressible  material.  Recent  numerical  experiments, 
reported  in  [6] ,  strongly  indicate  that  the  degeneracy  which  has 
led  to  the  use  of  such  techniques  as  reduced  integration  is  not 
present  in  the  case  of  the  p-version.  The  main  theorem  of  this 
paper,  Theorem  3.1,  verifies  that  conjecture. 

We  shall  briefly  illustrate  the  previous  discussion 
with  one  of  the  computations  from  [6] .  Figures  1  and  2  show 
the  relative  error  in  energy  versus  the  number  of  degrees  of 
freedom  (on  a  log-log  scale)  for  an  "edge  cracked  square  panel". 
The  mathematical  model  used  is  that  of  2-dimensional  plane 
strain  elasticity.  In  Fig.  1  the  Poisson  ratio  is  0.3  and  in 
Fig.  2  it  is  0.4999  (i.e.,  the  material  is  almost  incompressible). 
For  the  h-version  the  functions  are  piecewise  linear,  and 
the  p-version  is  based  on  piecewise  polynomials  on  the 
triangular  mesh  shown  in  the  upper  left  hand  corner.  By 
comparison  of  the  two  figures  it  is  evident  that  the  p-version 
is  quite  insensitive  to  changes  in  the  Poisson  ratio,  whereas 
the  h-version  for  v  =  0.4999  does  not  seem  convergent  at  all 
(to  be  more  precise,  one  is  not  in  the  asymptotic  range  even 
with  500  degrees  of  freedom) . 

A  remark  concerning  the  convergence  rates.  In  Fig.  1 
it  is  seen  that  the  p-version  has  a  convergence  rate  which 
is  twice  that  of  the  h-version.  This  is  a  consequence  of  the  fact 
that  polynomials  are  well-suited  to  approximate  a  wide  class  of 
singular  solutions  to  elliptic  boundary  value  problems 


Fig.  1.  Relative  error  in  strain  energy  vs.  number  of  degrees 
of  freedom. 

Poisson  ratio:  0.3. 

(Reproduced  by  permission  of  the  authors;  j.  Wiley  &  Son 
for  IJNME. ) 
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Pig.  2.  Relative  error  in  strain  energy  vs.  number  of  degrees 
of  freedom. 

Poisson  ratio:  0.4999. 

(Reproduced  by  permission  of  the  authors; 

J.  Wiley  &  Son  for  IJNME.) 
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as  the  degree  increases  (for  more  details  see  [7]  )  .  The 
convergence  rates  shown  in  Fig.  1  are  governed  by  the  singularity 
of  the  solution  at  the  crack  tip. 

The  analysis  presented  in  this  paper  does  not  study 
the  sensitivity  of  this  "corner  singularity"  convercence  rate 
to  changes  in  v.  Rather  the  results  are  directed  at  the  case 
when  the  convergence  rate  is  determined  by  singularities  of 
the  solution  inside  the  domain  (see  remarks  following 
Theorem  3.1). 

Numerical  evidence  (as  the  one  shown  here)  suggests  that 
the  "corner  singularity"  convergence  rate  is  also  insensitive 
to  changes  in  v  whenever  the  "corner  singularity"  itself  does 
not  change  with  v.  in  some  cases  however  (cf.  [6])  the 
"corner  singularity"  detrends  on  v,  and  therefore  naturally  also 
the  convergence  rate  of  the  p-version. 

The  contents  of  this  paper  are  as  follows.  In  section  2 
we  give  a  precise  statement  of  the  continuous  problem  to  be 
studied  in  this  paper  and  a  general  stability  analysis  for 
a  wide  range  of  approximate  solutions.  For  V  in  a  large  class 
of  closed  subspaces  of  H*(P)  we  assign  a  numbar  C  (  1/  )  that 
measures  the  stability  in  a  sense  made  precise  by  Lemma  2.1. 

In  section  3  we  establish  uniformly  valid,  optimal  rate  error 
estimates  for  the  p-version.  The  proof  of  the  main  theorem, 
Theorem  3.1,  relies  on  a  result  concerning  the  invertibility 
of  the  divergence  operator  as  established  in  [16] .  In  the 
appendix  we  prove  a  uniformly  valid  a  priori  estimate  for  the 
continuous  problem  which  is  needed  for  the  proof  of  Theorem  3.1. 


An  analysis  similar  to  that  presented  in  this  paper 
could  also  be  carried  out  for  the  Timoshenko  plate  model. 

Throughout  this  paper  H^(n)  denotes  the  standard  Sobolev 

2 

space  of  functions  that  have  all  derivatives  of  order  <  k  in  L  (ft) 
Whenever  no  ambiguity  is  possible  the  same  notation  shall 
also  be  used  for  vector  valued  functions  whose  components  are 
in  denotes  the  dual  of  H*(S2)  (and  not  of  Hk(fi)  as 

sometimes  used) .  Spaces  with  noninteger  exponents  are  defined 
by  complex  interpolation  (cf.  [9,13]).  The  norm  on  the  space 

g 

H  (Q)  is  denoted  I  •  II  n.  The  same  conventions  apply  tc  Sobolev 
spaces  on  the  smooth  boundary  3ft . 
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2 .  Formulation  of  the  continuous  problem  and  a  general 

stability  analysis. 

The  problem  which  we  shall  study  is  that  of  plane  strain  linear 

,  ,  .2  00 

elasticity  on  a  bounded  domain  ft  c  DR  with  a  C  boundary 

9ft-  Formulated  in  terms  of  displacements  uv =  (uv  ^,uv  ^ ) 

the  continuous  problem  is 

(1)  -  Auv  -  ~ V(V«uv)  =  F  in  n, 

where  F  represents  an  external  force  (to  be  specific  F  is  the 
external  force  scaled  by  2(1+v)/e).  0  <  v  <  i  denotes  the 
Poisson  ratio  and  0  <  E  the  Young  modulus;  a  value  of  v  close 
to  corresponds  to  a  nearly  incompressible  material.  For 
boundary  conditions  we  consider  the  natural  (stress)  boundary 
conditions  associated  with  (1) ;  these  may  be  written  as 
2  v 

<2)  ^£ij(-v)  +  6ij  TTfv  7*!2v)nj  =  gi  '  1  <  i  <  2  on  3S), 

2  2 
where  {n.}.  .  is  the  unit  outward  normal  to  3ft,  {£. . }r  .  , 

D  1=1  iD  i,D=l 

is  the  strain  tensor,  given  by 


e . 


ID 


-  ?  ( 


5S7 


,  d  a 

uj  +  HT  U1  > 


and  {g^}.  .  represents  the  force  being  applied  at  3f)  (here 
scaled  by  (1+  V)/E) .  The  boundary  value  problem  (1)  and  (2) 
has  a  solution  iff 


(3) 


[  F  R  dx  +  2  f 

V  3fi 


g  R  ds 


0 
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for  any  rigid  motion  R. 

In  the  following  it  shall  always  be  implicitly  assumed 
that  F  and  g  satisfy  this  relation.  The  solution ,  of  course,  is 
only  determined  up  to  a  rigid  motion. 

Let  denote  (H3-  (ft)  )  /{ rigid  motions)  with  the 

standard  norm  for  a  quotient  space  (also  denoted  "  * " i  * 
Solving  the  boundary  value  problem  (1)  and  (2)  is  equivalent 
to  minimizing  the  energy  expression 


(4) 


II 


l  (^(v)) 

Li'j-1  13 


1-2V 


(V-v) 


dx 


F  v  dx 


2 


j  g  v  ds 
dn 


in  (H^tfi))  ^  (or  to  make  the  solution  unique). 

Let  K  be  a  closed  subspace  of  tf3(ft)  (not  necessarily  proper) 
and  let  uv  ^  denote  the  minimizer  for  (4)  in  V .  For  V  =  H  iQ) 
we  drop  the  subscript  V  to  be  consistent  with  earlier  notation. 

The  specific  application  we  have  in  mind  regarding  numerical 
procedures  is  of  course  to  let  V  be  some  space  of  (relatively 
high  degree)  piecewise  polynomials  on  some  triangulation  of  Q. 

In  this  section  we  establish  an  estimate  on  uv~ uv  ^ 
which  is  valid  uniformly  in  v  (the  nontrivial  and.  interesting 
case  being  when  v  approaches  --) .  For  reasons  of  exposition 
we  have  chosen  to  give  this  result  in  a  general  form  here  and 
then  devote  the  next  section  to  the  application  for  the  p-version 
of  the  finite  element  method. 
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If  we  introduce  the  space  V  =  =  {V*v  |  v  e  l/}  and 

Lagrangian  multipliers  then  it  is  quite  clear  that  uv  y  may 
also  be  characterized  as  the  first  component  of  the  vector 
(v,r,s)  e  V  *  V  *  V  which  satisfies 


Let  C(l')  denote  the  infimum  of  all  constants  C  that  can  be  used 
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with  the  construction  (7) .  As  a  tool  for  our  later  analysis 
we  prove  the  following  stability  estimate  for  the  form  By . 

Lemma  2.1.  Let  By  denote  the  bilinear  form  defined 
in  (6).  Let  V  be  a  closed  subspace  of  with  the  property 

that  V  =  V»V  jt  0  is  closed  in  L^(fi). 

There  exists  a  constant  d  >  0,  independent  of 
v  €  ] 0,1/2 [  and  V,  such  that 

inf  sup  I  Bv  ( (v,r,s)  ;  )  |  _>  d  C(  V  )~2 

(v,r,s)eS  (w,<J> S 

where  5  =  {(v,r,s)  6  V*  V  x  V  |  lyl +  Irl Q  n  +  B  sB  Q  Q  =  1}. 

Proof :  One  argues  as  in  [4] ;  we  give  just  a  sketch  of 
a  proof  here  in  order  to  ensure  that  the  statement  about  the 
specific  form  of  the  lower  bound  is  true. 

Prom  Korn's  inequality  (cf.  [10]) 

(8)  E^((y,r,s);  (y,r,-s))  >  d(#  X* +  <  r*  jj  ,  . 

d  in  the  following  takes  the  role  of  a  generic  (small)  constant 
that  does  not  depend  on  v  or  V  .  If  w  is  selected  such  that 
7*w  =  -  s  and  flw8^  ^  <_  2C  (  V  )  I!  sO  q  ^  (this  can  be  accomplished 
according  to  the  definition  of  C(  1/ ) )  then 

(9)  Bv((v,r,s);  <w,  0 , 0) )  >  II  sll  2  ^  -  d_1C  (  V  )  II  yll  x  ,  QH  8|  Q  ,  Q 


12 


Taking  an  appropriate  linear  combination  of  (v,r,-s)  and 

(w ,0,0)  we  obtain,  due  to  (8)  and  (9),  a  (w,4>,<|')  that  satisfies 
the  estimates 

Bv((v,r,s);  (w,4>,*))  >  d(lvljfft  +  +  C  (  V)  "2I  si 


l-1  i,a  +  ,<*>,o,n  +  * ^ * o , n—  ,-,i,ft  +  Brl,o,ft  +  (1+c(  v  )  )Bsllo,ft 

It  is  easy  to  see  that  C  (  (/  )  is  bounded  away  from  0 
independently  of  1/  ,  and  that  the  estimates  (10)  therefore 
lead  to  the  desired  result.  0 


Based  on  Lemma  2.1  we  now  prove  an  estimate 
fOT  uv  -  uv  # V  • 

Lemma  2.2.  Let  uv  be  the  minimizer  of  (4)  in  W^fft) 
and  uv  y  the  corresponding  minimizer  in  the  closed  subspace  1/ 
Assume  that  the  space  P  =  V*  V  /  0  is  closed  in  L  (ft)  and  let 
C (  V)  denote  the  constant  defined  previously.  There  exists 
a  constant  D,  independent  of  v  €  ] 0 , 1/2 [  and  V  ,  such  that 

*-v  “  Sv,  V  *  i,n  +  i-2'v  11  V  *  ^-v  “  ^v,  P  ^ 11  0,ft 

<DC(1M2  (lnf^luv- vlliS1  +  inf  'j^r  7'Sv  -  0,n. 

Proof :  This  estimate  is  trivial  as  long  as  v  does  not 
get  close  to  1/2.  It  is  therefore  sufficient  to  consider 
1/4  <  v  <  1/2.  The  symmetric  bilinear  form  Bv  introduced 
before  satisfies 

(11)  |Bv((v,r,s)  ?  (w,<p,iM  )  |  I  d(iiYiii/q+ Br*0,ft+ ,s!!o,ft^ 

*  Csh.B  +  '♦'o,n  + 
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for  v  in  this  interval.  It  is  not  difficult  to  see  (cf.  [3]) 
that  Lemma  2.1  combined  with  (11)  leads  to  the  inequality 


2v  i/2  2v 

"Hv  ”  -v,^*  i,^  +  ^1-2V^  * ”  Hv  ,1^  B  o,n  +  IT’Tv  11  7‘  ^v~-v  ,  B  o,ft 


1.  (l+Dd^Cd/ftlinf^Hu^vlj^^  +  inf^  J  V2  V-u^-wI^ 

+  ^£fp  11 1-2V  V*-v  "  w*  0,n)  ' 

from  which  the  desired  estimate  follows  immediately.  □ 
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3.  The  p-version  F.E.M.:  uniform,  optimal  error  estimates. 

In  this  section  we  study  the  convergence  properties  of 
approximate  solutions  to  (1)  and  (2)  that  are  obtained  by 
minimizing  the  energy  expression  (4)  over  a  family  of  finite 
element  subspaces  associated  with  the  p-version. 

Let  1  be  a  triangulation  of  the  domain  Q . 

We  assume  that  no  vertex  of  a  triangle  lies  in  the  interior  of 
a  side  of  another.  Triangles  of  1  may  possibly  have  one  curved  side 
(in  common  with  3ft) .  We  assume  that  any  of  the  triangles  with  one 
curved  side  can  be  mapped  affinely  onto  a  Lipschitz  domain  of 
the  form 

{(xlfx2)  |  0  <  x1  <  1  ,  0<x2<  f  (x1)} 

where  f  is  monotone  and  satisfies  f (0)  =  1  ,  f(l)  =  0  . 

Let  denote  the  set  of  continuous  functions  whose 

restriction  to  each  triangle  is  a  polynomial  of  degree  at  most  p. 
pip]/  1  denotes  the  set  of  functions,  whose  restriction  to  each 
triangle  is  given  by  a  polynomial  of  degree  at  most  p  (no  inter¬ 
element  continuity  requirements),  with  the  additional  property  that: 


R.l.  "At  any  internal  vertex  (x®,x2)  of  the  triangula¬ 
tion,  where  four  triangles  meet  as  the  intersection 
of  two  straight  lines,  the  corresponding  values 
(x^,x2) }f_i  satisfy 


l 

i  even 


l 

i  odd 


Qi(x? 


If 


(cf .  Fig.  1  in  [16] ) . 

The  space  (P ' °) 2  /  {rigid  motions)  is  also  for  short 
denoted  (P^'°)2,  p  >  i. 
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For  the  proof  of  the  main  theorem  in  this  section  we 
need  the  following  3  auxiliary  results. 

Lemma  3.1.  Let  c((P^'^)2)  denote  the  constant 
introduced  in  section  ?..  There  exist  constants  C  and  K  , 
independent  of  p  ,such  that 

c(  (P[p]  »0)2)  <  c  pK 

for  any  p  :>  1. 

Proof:  This  follows  directly  from  Theorem  2.1  and 

Remark  2.2  in  [16]. 


Lemma  3.2.  Let  k  >  1.  There  exists  a  constant  C, 


such  that 


inf  !v  -  QPH.  <  c.  p“k+1ly«,  a 

qP6  (p[p],0}2  1,n  k  k.n 


k  2 

for  any  p  >  1  and  any  v  6  (H  (fi))  • 


Proof :  See  [7] . 


Lemma  3.3.  Let  k  ^  0.  There  exists  a  constant  such 
that  for  any  p  :>  1  and  any  w  €  Hk(fl)  , 


inf  Iw  -  QPl0((!  <  Ckp-klwlk(!i  , 

Qp 

where  the  infimum  is  taken  over  qP  e  V»((P^'®)^). 

Proof :  In  [16]  it  is  shown  that  V*((P^'®)2)  =  ptp-ll  ,-l 

for  p  sufficiently  large.  Due  to  the  inclusion 

p[p  1],0  ^  p[p  1],  1  tjje  estimate  follows  immediately  from  [7] 


(as  did  Lemma  3.2). 
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We  are  now  in  a  position  to  prove  the  main  result  of 
this  paper,  a  theorem  that  provides  optimal  error  estimates 
for  the  p-version  of  the  finite  element  method  based  on  direct 
minimization  of  the  energy  in  displacement  form  (4) .  The 
important  feature,  that  distinguishes  these  estimates  from 
the  optimal  ones  that  hold  for  the  more  standard 

h-version,  is  that  they  are  uniform  in  v.  We  shall  give  a 
discussion  to  clarify  in  what  sense  these  estimates  are  optimal 
immediately  following  the  proof. 


Theorem  3.1.  Assume  that  s  >_  -1  and  v  e  ]0,l/2[. 

Let  uv  denote  the  solution 

to  (1)  and  (2)  .  For  any  p  >_  1  let 

uv  p  denote  the  element  (or  rather  the  equivalence  class) 
in  that  minimizes  (4).  Given  e  >  0  there  exists 

Cs  P  (depending  on  s,  e  but  independent  of  v  and  p)  such 
that  whenever  F  €  HS(ft)  and  g  €  Hs+^//^(3^) 


*Hv  "  2?.vfp*i,n  +  (i-2v)  ^  1  v*  ^-v  “  1  o,n 


iCS,e  P 


- s-l+e 


('£'  s,n  +  "2s s+i/2,aoJ 


Proof:  By  Lemmas  2.2  and  3.1 

“  -v,p*  l,fl  +  1-2U  1  v*  ('iv  ~  11  o,n 

(15)  > 

<Cp2K  inf  luv-Yllfn  +  inf  1  r=2v  ' 

'  v  w 

where  the  first  infimum  is  taken  over  v  e  (p(P)'0)2  an<3  the 
second  over  we  V •  ( ( P • 


An  application  of  Lemmas  3.2  and  3.3  show  that  the  right 
hand  side  is  bounded  by 


CN  ? 


2K-N- 


(S-VllN+2,f2  +  1-2V  11  V*HvBN+l,n,  ' 


and  because  of  the  uniform  a  priori  estimates  given  in 
Theorem  A.l  of  the  appendix,  this  is  at  most 


CN  P 


2K-N- 


N_1  +  *«* 


N+l/2,3nJ  * 


If -we  replace  1/(1-2V)  in  the  left  hand  side  of  (15)  by 
1/2 

(1/(1”2V))  (which  is  certainly  smaller)  we  conclude  that 

i  1/2 

'Sv-Sv.p'u*  (i=jv>  ,7-<Hv-av,p>'o,n 


iCNP 


2K-N- 


‘(•'Vo  *  'sVi/2 


,3!i)  • 


Since  Hv,p  can  also  be  interpreted  as  the  projection  of  uv  onto 
t  p  1  0  2 

(p  '  )  in  the  energy  inner-product  associated  with  (4) , 

i  1/2 

'-v  "  ^v,P' i,n  +  11  v*  “  Uv,p) ffo,f2 

(17)  <  c(iav«lr0  ♦  (r^v)1/2  »7.uv»0ffi) 


(the  last  inequality  is  due  to  Theorem  A.l).  Interpolating 
between  (16)  and  (17)  subject  to  the  finite  dimensional  constraint 
that 
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ll  l  E 


R  dx  dy  +  2  g  R  ds  =  0 


{  2  R 


for  any  rigid  motion  R  (cf.  Theorems  12.5  and  13.3  of  [13]) 


we  get  that 


1  -v  "  -v^1  l,n  +  ^l-2v^ 


6(2K-N-1), 

iCN,ep  I'V 


(N+l)0-l/n+  '211  (N+l)6  -1/2, 39. 


Choose  N  larger  than  2K(s+l)/e  and  define  9  to  be  (s+1)/(n+1). 
This  ensures  that  2K9  <  e,  and  from  the  inequality  (18)  we 


thus  get 


1 

lcs,eP's'1+cCE's,n  *  •j'.+i/j.aa)  ■ 


exactly  as  desired. 


Remark  3.1.  The  estimate (14)  is  optimal  in  the  following 

(P-i  1  t  0  9 

sense:  If  there  exists  vQ  and  a  sequence  6  (p  J  )  , 

pj+l^Pj  —  ^  <  00  *  Pj  +  00  »  such  that 


-s-1-6 


'Sv„  -  Xj'i.a  <  tpj 


for  some  <5  >  0,  then  one  may  conclude  (cf.  [7])  that 


u„  e  H?+2+t(T.) 

loc  x' 
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for  any  t  <  6  and  any  triangle  £  Z  (u^  €  ^(T^)  means 

that  uv  e  Hs+^+^(K)  on  all  compact  subsets  K  of  T\) .  This 
immediately  implies  that 


(19)  F  €  H^(T.) 

for  any  triangle  7\  e  E.  For  a  general  F  6  Hs(fl) ,  (19)  is 
obviously  not  satisfied  and  the  best  possible  convergence 
rate  is  therefore  p”S_^. 

For  a  particular  class  of  singularities  typically 
associated  with  corners  of  the  domain  ft  it  has  been  verified 
that  one  does  better  than  (14) ,  indeed  the  p-version  converges 
at  twice  the  r?  *  of  the  h-version  finite  element  method 
(cf.  [7]).  As  already  observed,  numerical  evidence  indicates 
that  this  improved  convergence  rate  is  also  insensitive  to 
changes  in  v.  This  claim  however  is  not  verified  by  the 
present  analysis. 


Remark  3.2.  The  norm  on  the  left  hand  side  of  (14)  is 
equivalent  to  the  energy  norm 


1 


( e  i  j  ( v ) )  2  + 


V 

1-2V 


(V*v) 


1/2 


with  constants  that  are  independent  of  v.  Let  Wv  denote  the 
energy  expression  given  in  (4) .  From  the  observation  about 
equivalence  of  norms  and  the  estimate  (14)  it  follows  that 
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o  i  V“v 


p)-W  I  Ck/L 


-2k-2+e 


The  situation  is  as  already  mentioned  quite  different  for 
the  h-version  of  the  finite  element  method.  If  the  minimization 
is  performed  over  all  continuous  piecewise  polynomials  of  degree 
<  k+1  on  a  triangulation  of  mesh  size  h,  then  we  have  the 
estimate 


(20)  0  <  Wv(u") -Wv(uv)  <Ck(Vh 


2k+2 


f'E'k.o  +  '2'iLi/j.an  ) 


for  any  fixed  value  of  vs  )0,l/2[.  The  estimate  (20) 
however  is  not  valid  with  a  constant  that  is  independent  of  V; 
rather  v  tends  to  00  as  v  approaches  1/2.  Numerical  evidence  (e.g. 
the  figures  shown  in  [15) )  suggests  that  the  best  possible  estimate  is 

0  <  Wv(uJ)-  Wv(uv)  <  Ckh2k 

with  independent  of  v.  That  is  to  say,  the  uniformly  valid 

convergence  rates  for  the  energy  are  reduced  by  a  factor  of 
-2 

h  from  the  optimal  ones.  Similar  statements  seem  to  hold  for 
the  convergence  rates  for  the  displacements. 
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Appendix.  Uniform  a  priori  estimates  for  the  continuous  problem 


Because  of  the  regularity  properties  of  solutions  to 

elliptic  boundary  value  problems  on  smooth  domains  it  is  clear 

k+2 

that  uv  ,  the  solution  to  (1)  and  (2),  is  an  element  of  H  (ft) 

provided  F  e  H^(ft)  and  g  s  H^+^//^(8ft).  In  this  appendix  we 

k+2 

shall  prove  that  the  norm  of  uv  in  H  (ft)  and  the  norm  of 
1  k+1 

J-J2V  v*Hv  11  (A)  stay  bounded  independent  of  v  (this  is 

exactly  the  result  which  was  used  in  the  proof  of  Theorem  3.1). 
This  result  may  seem  surprising  at  first,  since  (1)  and  (2)  looks 
like  a  singular  perturbation  problem  for  v  close  to  1/2;  what 
we  therefore  in  effect  show  is  that  there  are  no  boundary 
layers. 

Let  us  introduce  the  following  auxiliary  problem 


(21) 


V  AUV  +  VPv  =  F 

in  ft 

V*yv  =  G 

in  ft 

with  boundary  conditions 
2 

(22) 


J,  "j  '  pvni 

3=1  J  J 


g^  ,  1  5  i  £2  on  3ft. 


This  is  a  sort  of  generalized  Stokes  problem;  we  now  proceed 
to  establish  apriori  estimates  for  and  in  terms  of 
F,  G  and  g  .  In  doing  so  we  rely  on  the  paper  of  Agmon, 

Douglis  and  Nirenberg  [1] .  Our  proof  therefore  consists  merely 
of  checking  that  the  system  (21),  (22)  falls  into  the  general 
class  studied  in  that  paper. 
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Lemma  A.l.  Assume  that  1/4  <  v<  1/2  and  that  k  is  an 

integer  >  0.  If  Uv  e  h2(£),  Pv  e  H1  f )  denotes  a  solution 

~  y  -  Vii 

to  (21)  and  (22)  for  data  F  6  H  (B) ,  G  e  H  (0)  and 
g  E  Hk+1^2(30)  then  Uv  £  Hk+2(n),  Pv  £  nk+1(S!)  and 


l2v'k+2,S!  +  ,Pv'k+l,B  i  Ck('£'k,B+  '5,k+l.a  +  "i'k+l/Z.JB 


+  'Sv'o.B  +  lpv'o,fl)  ■ 


with  C.  independent  of  v . 


Proof :  If  we  denote  the  vector  (UV,PV) 

by  Vv  then  (21)  may  be  written 


~  -va 

v" 

F 

Vv  = 

V* 

0  _ 

—V 

_  G_ 

The  symbol  of  this  system  is 

*(£) 

By  defining  t^  =  2,  j  *  1,2,  t3  =  1 , 
and  s3  =  -1  we  obtain 


o 

0  -v|5|2  52 

?  i  5  0 


deg  (J-^j  (£) )  <  +  t .  . 


s . 
i 


0,  i 


1,2, 


Since  in  this  proof  we  follow  the  notation  of  fl]  exactly,  we  shall 
not  care  to  explain  the  meaning  of  these  indices  but  rather  refer  the 
reader  to  that  paper  for  more  details.  The  principal  part  8.!.  (£)  is 
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equal  to  and  we  compute 

L(§)  =  det  (*[..(£))  -  vl£|4  • 


which  ensures  the  uniform  ellipticity  of  the  system  (condi¬ 
tion  1.7  on  page  39  of  [1]). 

If  £  and  i *  are  two  nonzero  real  vectors  then  the 

equation 


(23) 


L(£  +  T£')  -  0 


has  the  roots  _ 

*  -  e-L'  1  i  /|5|2|5'|2-  (£•£■>* 


both  with  multiplicity  2.  If  £  and  £'  are  linearly  independent 
then  |£|2|£' 1 2  “  (£*£') 2  >  °»  and  we  conclude  that  the 
equation  (23)  in  this  case  has  exactly  2  roots  (counted 
according  to  multiplicity)  with  positive  imaginary  part. 

This  verifies  the  supplementary  condition  on  page  39  of  (11. 

The  boundary  condition  ( 22)  may  be  represented  by  the 
symbol 


^l^l"*"  2  ^2n2  2  ^ln2  ~ nl 

*2  ^2nl  ^2n2+  "2  ^1^1  “n2 

(The  unit  normal  n  depends  on  Q,  the  particular  point  at  the 
boundary  of  0.) 

By  choosing  r.  =  -1,  h  *  1,2,  we  obtain  the  inequalities 


§(  Q,£)  = 


deg  (Bhj  (  0,£))  <  rh  +  t..  , 
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and  since  these  actually  are  identities  we  conclude  that 
the  principal  part  B*  is  identical  to  B.  The  product 
B  (  Q  ,C)  *4  (£)  is 

“  vlil  ^lnl+2  “  2  I— I  ^2”  ^"l  ^lnl+^1^21^ 

--  2  lii%VSln2  -»ltl2(52r|-«lnl)-i:2n2 

For  the  complementing  boundary  condition  (p.  42  of  [1]) 
to  be  satisfied  we  have  to  check  that  every  nontrivial  linear 
combination  of  the  two  rows  in  the  matrix 

B  (  Q  ,  t+ t  n)  •  1  ( fc+  t  n) 

is  nonzero  modulo  the  polynomial 

M+ (  Q , t, t)  =  (t  •  Tq (t»n) ) ^  . 

This  has  to  be  satisfied  for  any  tangent  vector  t  /  0. 

If  we  consider  *  (first  row)  +  C2  x  (last  row)  the 
third  component  has  the  form 

(C*  t  +  t  C  •  n)  r 

where  C  is  the  vector  (C^,C2).  This  can  be  0 

modulo  M+( Q  ,  t, t)  only  if  C»t  =  C*n  =  0,  i.e. , iff  C  =  0. 

We  have  therefore  verified  the  complementing  boundary 
condition. 

Since  \>  is  always  in  the  interval  ] 1/4, 1/2 (  (i.e.,  the 
lower  bound  in  the  uniform  ellipticity  estimate  does  not 
degenerate)  the  desired  estimate  follows  directly 
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from  Theorem  10.5  of  [1].  □ 

Based  on  Lemmas  2.1  and  A. 1  we  now  prove  the  following 
apriori  estimates  for  uy  ,  the  solution  to  (1)  and  (2). 

The  important  feature  is  that  the  constants  are 

independent  of  v  £  ]0,l/2[. 

Theorem  A.  i.  Let  k  be  an  integer  >_  -1. 

If  uv  £  H^(n)  denotes  the  solution  to  (1)  and  (2) 

for  data  F  £  Hk(n),  g^€  Hk+1//2On)  then 

uv  £  Hk+2(n)  and 

(24)  lnvlk+2,n  +  1^27  * 7 v* k+i , ft  -  Ck(,£,k.n+  S2-l!k+l/2, an) 
with  independent  of  v. 

Note;  The  solution  to  (1)  and  (2)  is  only  uniquely 
determined  up  to  a  rigid  motion  (i.e.,  in  H*(n)).  The 
first  norm  on  the  left  side  in  (24)  refers  to  the  norm 
of  the  corresponding  equivalence  class. 

k+?  k 

Proof :  That  £  H*  'in)  provided  F  £  H  (Q)  and 

2_  e  Hk+1//2On)  follows  directly  from  standard  regularity 
results  for  elliptic  boundary  value  problems  on  smooth 
domains.  That  the  estimate  (24)  is  uniform  in  v  for 
v  £  ]0,c[,  c  <  1/2,  is  also  clear  from  that  theory.  We  shall 
now  show  that  (24)  is  valid  uniformly  in  v  e  ]c,l/2[, 

for  an  appropriate  choice  of  c»l/4. 

In  the  case  k  -  -1  this  result  follows  immediately 

V  =  To  be  more  specific 


from  Lemma  2.1,  by  taking 


r 

^  6 


I  u  8 


2V 


Vl,«  +  5757  "  V-Hv“  0,n 


<  d_1C(H1)2  sup  |B  ((u  ,r,s);  (w,4>,4>))! 

(w,4>^)e  s 


with  r  =  1//2  7*uv  and  s  =  V »uy .  The  right  hand 

side  is  according  to  (5)  equal  to 


d_1C ( H  X)2 


|  F  w  dx  +  2  j  £  w  ds|  , 


+  2 

n 

and  this  gives  the  desired  estimate. 

In  the  case  k  >  0  let  us  introduce 


U  a  U 
—v  — \> 


Pv  "  l-2v  V'-v  ' 

then  Uv  €  H2(n),  P^  e  H1(n)  satisfies  (21)  and  (22)  with 
F  =  v  F  ,  G  =  V *u  and  2.  “  2.*  By  Lenuna  A.l, 


Ju  I 


B  7  •  u  B 


-V  k+2,n  i-2v  -v  k+l,n 


(25)  =  BU  |  +  I P  B 

—v  k+2,J2  v  k+l,fl 


<  Ck(v,F|k^n+  ,7*)iv)l'k+lfn+  B  2®  k+1/2/  3  0+  l-2v  ^'Hv*  0,^* 


The  last  two  terms  on  the  right  hand  side  can  by  an  applica¬ 
tion  of  this  theorem  with  k  = -1  (which  has  already  been  proven) 


be  estimated  by 


Ck^l-l-i,fl  +  ,2' -1/2, 3^  ' 

and  therefore 

'-v'k+a,^  +  I^2v  lv*iJv,lk+i,n 

ick(»Flk>fi  +  *v*uv*k+lfn  +  *2lk+i/2,3n^  • 

Since  v  e  ]c,l/2[,  we  may  by  choosing  c  sufficiently 

1  v 

close  to  1/2  arrange  that  <  2  *  i-2v  '  anc*  t^us 

absorb  the  term  involving  V*uv  on  the  right  side 
into  the  corresponding  term  on  the  left.  This  concludes 
the  proof  of  the  theorem.  □ 
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